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Analyzing real data on international trade covering the time interval 1950-2000, we show that 
in each year over the analyzed period the network is a typical representative of the ensemble of 
maximally random weighted networks, whose directed connections (bilateral trade volumes) are 
only characterized by the product of the trading countries' GDPs. It means that time evolution 
of this network may be considered as a continuous sequence of equilibrium states, i.e. quasi-static 
process. This, in turn, allows one to apply the linear response theory to make (and also verify) 
simple predictions about the network. In particular, we show that bilateral trade fulfills fluctuation- 
response theorem, which states that the average relative change in import (export) between two 
countries is a sum of relative changes in their GDPs. Yearly changes in trade volumes prove that 
the theorem is valid. 
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I. INTRODUCTION 

In recent years, an extensive research effort has been 
devoted to analyzing the structure, function, and dvnam- 
ics of the international trade network (ITN) [iH^ from 
a complex network perspective The knowledge of 
topological properties of this network and its evolution 
over time is not only important per se (e.g., because it 
enhances our descriptive knowledge of the styhzed facts 
pertaining to the ITN), but it may also be relevant to a 
better explanation of macroeconomic dynamics (Tol4T7|. 

Here, we use quantitative and numerical (data-driven) 
methods originating from statistical mechanics to de- 
scribe and predict the behavior of ITN. We analyze a 
set of year-by-year trade relationships between all coun- 
tries of the world, covering the time interval 1950 — 2000. 
Although the total number of countries and the over- 
all economic conditions influencing the network change 
over the course of the period, in each year ITN is shown 
to be a typical representative of the ensemble in which 
every network, is assigned the probabihty [3, 
P{G) oc e-^(^); where H{G) = plays the 

role of network Hamiltonian; Wij represents the volume 
of trade between two countries, i and j; Oij oc {xiXj)~^ is 
the field parameter conjugated to this trade connection; 
and XiXj corresponds to the product of the GDPs of the 
trade partners. 

Behind the descriptive power of our approach (which 
has been confirmed in a number of tests reported in this 
article consisting in comparison of GDP-driven Monte 
Carlo simulations of the trade network with real data on 
ITN), it is also important to stress the predictive abili- 
ties of the model. In particular, we show here that bilat- 
eral trade fulfills a simple fluctuation-response theorem 
p3 |. Supported by the well-known qualitative flndings 
about economic crises, we argue that the theorem pro- 
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vides valuable quantitative insights into dynamics of ITN 
in the time of crises, and, in the future, may be used to 
uncover mechanisms underlying the emergence of world- 
wide crises. 



II. DATA DESCRIPTION 

The results reported in this work are based on the trade 
data collected by K. S. Gleditsch [20] that contains, for 
each world country in the period 1950 — 2000, the de- 
tailed Hst of bilateral import and export volumes. The 
data are employed to build a sequence of matrices, W(t), 
corresponding to snapshots of weighted directed ITN in 
the consecutive years, t = 1950, .... In the network, 
each country is represented by a node and the direction 
of hnks follows that of wealth flow. The entry, Wij{t), of 
the trade matrix, W(t), represents the weight of the di- 
rected connection. From the point of view of the country 
denoted by z, Wij{t) refers to the volume of export to j, 
while, from the point of view of the country labelled by 
j, it is seen as the volume of import from i. Precisely due 
to differences in reporting procedures between countries, 
when analyzing the data one often encounters small de- 
viations between exports from i to j and imports from i 
to j. To overcome the problem, in our analysis we deflne 
Wij{t) as the arithmetic average of the two values. 

In this article, apart from trade matrices, which con- 
tain complete but often excessively detailed information 
about ITN, we also use several other quantities that make 
theoretical description of the network possible. In par- 
ticular, to characterize the economic performance of a 
country we use its total GDP value, Xi{t). To get the 
whole set of total GDPs, {xi{t)}, we simply multiply 
the GDP per capita by the population size of each coun- 
try [21]. Furthermore, to describe the intensity of the 
trade relationships of a country we deflne the so-called 
out-strength, s°'"*(t), and in-strength, s-'^(t), of the cor- 
responding node. The quantities are calculated as the 
total weight of connections (outgoing and incoming, re- 
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spectively) that are attached to the node and they repre- 
sent total volumes of export and import of the considered 
country in a given year, t. 

All the data used in this study are given in milhons 
of contemporary U.S. dollars. The disturbing effects of 
inflation are ruled out in a natural way by the fact that 
whenever the variables, Xi{t) or Wij{t), are used in the 
calculations, they are intrinsically divided by the normal- 
ization constant that equals the sum of all variables of a 
given type, i.e. X{t) = J^i^M or T{t) = Ei,jW^ijW) 
respectively. (In what follows, whenever there is no con- 
fusion we will often omit explicit time dependence of the 
quantities). As a byproduct, the above observation al- 
lows one to present the main results of this article in 
a very concise way, with the help of relative quantities 
defined as follows: 



this ensemble can be written as 
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III. EXPONENTIAL RANDOM GRAPHS: 
STRUCTURAL HAMILTONIAN OF THE 
INTERNATIONAL TRADE NETWORK 

The first contributions dealing with international trade 
from a complex network perspective used a binary- 
network approach, in which one has assumed that a (pos- 
sibly directed) link between any two countries is either 
present or not, depending on whether the trade volume 
that it carries is larger than a given threshold . With 
reference to this hue of research we would hke to highhght 
the paper by Garlaschelli and Loffredo (2004), published 
in Physical Review Letters Q- In the paper, the authors 
used the same real-world data to analyze an unweighted 
and undirected version of ITN, i.e. a network of partner- 
ship in trade. They have shown that the total GDP of a 
country, can be identified with the fitness variable [22| 
that, once a form of the probability of trade connection 
between two countries is introduced, completely deter- 
mines the expected structural properties of this network. 
This in turn implies that ITN viewed as a binary network 
is a typical representative of exponential random graphs 
[Tsl . Furthermore, one can expect that the same holds 
true for the weighted version of this network [l^, • 

To verify this conjecture, we start by considering the 
most general ensemble of directed weighted networks, 
which is described by the following Hamiltonian 



(2) 



with a separate parameter 9ij coupling to each weighted 
connection. Our aim is to check whether the Hamiltonian 
is correct and, if so, how the parameters {9ij} depend 
on different indicators characterizing the global economy. 
To do this, we first examine the ensemble as it stands. 

Thus, given that Wij is a real number greater than 
(as is true for trade volumes), the partition function of 



'dwii = \{\{j-. (3) 



This allows us to rewrite the probability of a network as 
P(G) = e-^(«)/Z = n.n,^,K., where 



Pij(wij) = e~ 



(4) 



is the probability that there is a directed link of weight 
Wij from i to j. The expression for pij that we arrive at 
is the exponential distribution. Its mean value. 



(5) 



can be used to calculate the average values of a node's 
out- and in-strength 



At this stage one can start to make comparisons of 
theoretical predictions with the empirical data on inter- 
national trade. With good reason, it is convenient to 
begin by putting together Eq. (O and the corresponding 
empirical relations (see Fig. [T]a, b, c, and d): 



(sr*)=^a;, and {s^) = Axi, 



(7) 



where A is the time-dependent parameter having the 
same value for both out- and in-strength of the nodes. 
Analyzing the expressions, one finds that the simplest 
way to merge the theoretical approach with real data is 
to assume a multiplicative form of the parameter 9ij^ i.e. 

Oij = 9i6j, (8) 

where 6i and 6j represent some single-node characteris- 
tics controlhng for the potential ability of the two nodes 
to be connected. One should note that the symmetric ex- 
pression for Oij, Eq. JH]), is consistent with observations 
made by other authors, showing the symmetric character 
of bilateral trade relations (see e.g. [3, 8]). 

To calculate explicit values of all the parameters {^i}, 
one just has to insert Eq. JH]) into the theoretical for- 
mula for (s°"*), Eq. ([6]), and then equate the obtained 
relation to the empirical one, Eq. (|7|). (The analogous 
calculations can be done for (s^").) As a result, one 
gets the following expression: (s°"*) = Oj^ /6i = Axi, 
which, when summed over z, yields an important relation 
between theoretical and empirical quantities describing 
ITN: T = = AX, from which it follows that 
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where ^i has been introduced earlier, in Eq. ([T}. Now, one 
can insert (HJ into (|5]). As a result one gets a stochastic 
version of the gravity model of trade |25l427| , 



T 



(10) 
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In comparison to the standard gravity equation for bilat- 
eral trade volume between i and j, 

= """'^m^y ^^^^ 

where f{Dij) represents transport costs in international 
trade, which depend on geographic distance between the 
two trading countries. 

The expressions, Eqs. (O, together with other relative 
parameters defined in Eq. ([T]), can be used to rewrite the 
most important results of this approach in a very concise 
way. In particular, as described in terms of trade, the 
average out- and in-strength of a node, Eq. ([6]), when 
divided by the world's trade volume, T, turns out to be 
equal to the country's share in the world's GDP, X, i.e. 

(a™*) = wr) = (12) 

In a similar fashion, the average weight of a directed con- 
nection when divided by T, is given by 

{V^J)=^^^J=^/V^J■ (13) 

Finally, the structural network Hamiltonian, Eq. ([2]), 
when written in relative variables has the following form 

HiG) = Y^Y.v^jV.j. (14) 

i j^i 

To verify the correctness of the assumed network 
Hamiltonian, a series of data-driven Monte Carlo sim- 
ulations employing the Metropohs algorithm has been 
performed. The obtained results, for four years: 1973, 
1974, 1975, and 1995, are shown in Figure [T] (diagrams 
placed in the same row refer to one year marked in the left 
chart). The first three years correspond to the 1973 oil 
crisis which caused high inflation and a global recession 
which affected all aspects of living within the 1970s. The 
figure shows that global economic crises, such as the one 
triggered by the 1973 oil crisis, do not affect the quahty 
of our approach. In particular, in panels e, f, g, and h, in 
the Figure [T] sets of all bilateral trade volumes recorded 
in a given year versus the product of the trading coun- 
tries' GDPs is compared with the corresponding set of 
weights of directed connections in a typical network of 
the considered ensemble. Although the two sets (clouds) 
of points are quite dispersed, they overlap significantly, 
and their shape is well-described by Eq. ([T3|) . Moreover, 
as shown in panels i, j, k, and 1 in the same figure, the 
distributions of trade volumes within these clouds fit very 
well with each other and agree with the distribution of 
expected trade flows, P{{vij)), testifying in favor of our 
simple, yet realistic, approach. 

IV. SUSCEPTIBILITY OF TRADE VOLUME TO 
CHANGES IN GDP: FLUCTUATION-RESPONSE 
RELATION FOR ITN 

Extensive comparisons between real data on interna- 
tional trade and its GDP-driven Monte Carlo simula- 
tions show that, although the total number of world's 



countries, N{t), and their GDPs, {xi(t)}, change over 
the analyzed period of 50 years, ITN is continuously 
well-characterized by the same Hamiltonian. This means 
that the time evolution of this network may be consid- 
ered as a continuous sequence of equihbrium states (i.e. 
quasi-static process) that is yearly sampled by the na- 
tional reporting procedures. Furthermore, since differ- 
ences between snapshots of ITN in the consecutive years 
are rather small, one can expect that they could be de- 
scribed with the help of liner response theory, of which 
the simplest (but not yet trivial) result is the fluctuation- 
response theorem [24]. 

In the case of exponential random graphs with the 
Hamiltonian given by Eq. (|14p the fluctuation-response 
theorem has the following form 

(vfj) - {vtjf = -d{v,j)/dr],j. (15) 

The l.h.s. of this expression describes fluctuations in rela- 
tive weight, Vij = Wij/T, of the directed connection from 
i to j, whereas its r.h.s. characterizes susceptibility of 
this hnk to its conjugated local field, ijij. The suscepti- 
bility is defined as the derivative of (vij) with respect to 
r]ij and describes what happens with the expected trade 
volume, (vij), when one changes the parameter rjij, which 
determines external conditions for the bilateral exchange. 

Taking into account the exponential distribution of 
weights, p{vij) = e~'^'^'"'^r]ij, which follows from Eq. (|4|) 
and whose variance equals the square of the mean, the 
fluctuation-response theorem, Eq. ([15]), can be trans- 
formed into the formula 

{vij) rjij 

Written in such a way, the theorem states that rela- 
tive changes in normalized (i.e. divided by T) bilateral 
trade volumes can be estimated on the basis of changes 
in the GDP of trade partners. Yearly changes in im- 
port/export volumes between different countries prove 
that the fluctuation-response theorem for ITN is correct 
(see Figure [2]). Although, we have obtained Eq. ([T6)) as 
the fluctuation-response relation for exponential random 
graphs, it is obvious that it can be also obtained as the 
logarithmic derivative of Eq. ([T3)) . 

Relying on Eq. ([T6|) one can, for example, expect that a 
decline of, say, 2 percent in the relative GDP of a country, 
given that its trade partners do not change their share 
of the world's GDP, will translate into a similar decline 
in all its bilateral trade volumes. The example shows 
that the theorem can be used to make simple predictions 
about the world-wide diffusion of trade-based economic 
perturbations. Furthermore, assuming that the struc- 
ture of ITN is a proxy for meaningful financial linkages 
between countries, the expression may also help one to 
understand how financial ripples originating in one coun- 
try propagate to other countries, giving rise (or not) to 
global financial contagion. In particular, in the light of 
this theorem, it becomes apparent that a crisis is ampli- 
fied if the epicenter country is better integrated into the 
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trade network [Isl. This happens because the dedine of 
its GDP, through the proportional downward effect on 
trade and financial linkages, affects more countries. It 
also becomes clear that the impact of a crisis on any 
target country is cushioned if the country in question is 
better integrated into ITN, due to the fact that, in the 
case of such countries, the dechne of only one bilateral 
trade volume does not significantly influence GDP. 

V. CONCLUDING REMARKS 

The current economic crisis illustrates a critical need 
for new and fundamental understandings of the struc- 
ture and dynamics of economic networks. This sentence 
opens the perspective article from the special issue of Sci- 
ence [9] on complex systems and networks entitled " Eco- 
nomic networks: The new challenges" by F. Schweitzer et 
al. In their article [llj, the authors summarize what we 
know and what we need to know about different economic 
networks (including ITN) to reduce the risk of global 
depression and to design effective strategies to promote 



economic recovery [28|. Our approach to ITN is in hue 
with this challenging research area. Having the mathe- 
matically tractable yet realistic model of ITN introduced 
here, and given its quasi-static time evolution, we be- 
lieve that, apart from the fluctuation-response theorem, 
Eq. ()16p, other well-known results of non-equilibrium sta- 
tistical physics [29] may be apphed to estimate recession 
(or economic growth) impact on international trade. 

The exponential random graphs are in common use 
within the statistics and social network analysis commu- 
nities as a practical tool for modeling networks for a few 
decades. Although such modehng gives a good quahta- 
tive description of the data, to our best knowledge, there 
is no evidence so far that with such approach one could 
achieve excellent quantitative accuracy and reproducibil- 
ity of real world phenomena. 
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FIG. 1: Structural properties of single snapshots of ITN. Diagrams placed in the same row refer to the same year. The 
data are shown in two ways, by using regular and relative quantities, cf. Eq. ((TJ). a, b, c, d, Total import and export volumes 
of all world countries in 1973, 1974, 1975, and 1995 vs. GDP (filled and open points correspond to imports and exports, 
respectively) and their comparison with the expected values described by Eqs. ((7)) and (p^ (sohd straight lines), e, f, g, h, 
Bilateral trade flows in the following years vs. the product of the trading countries' GDPs (points) as compared with their 
theoretical prediction based on Eqs. (jS)) and (|13|) (line). Black points correspond to real data, while gray points represent trade 
volumes obtained from GDP-driven Monte Carlo simulations. In the insets of all the panels (e, f, g, and h), comparison of the 
expected theoretical values (straight lines) to mean values of the real data is shown. Since trade flows smaller than a given 
threshold are rarely specified in economic reports (in particular, the considered data set ^] does not contain trade volumes 
smaller than 1000 USD), the cloud of black points cover smaller area than the one corresponding to numerical simulations, i, 
j, k, 1 Distributions of trade volumes in the considered years. Filled and open squares correspond to real and simulated data, 
respectively. The solid lines represent distributions of expected trade flows which, for each pair of countries in a given year, 
can be calculated using Eqs. (jS)) or (|13|). 
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FIG. 2: Fluctuation-response theorem for ITN. To confirm validity of the theorem, for each connection in the trade 
network at the turns of 1950 — 2000, one has calculated quantities corresponding to both sides of Eq. (|16() : the relative change 
in normalized trade volume, dvij/vij c:^ {vij{t-\- 1) — Vij{t)) /vij{t), and the sum of relative changes in GDP of trading countries, 
d^i/Ci + d^j/Cj^ where d^i/^i c=i {^i{t -\- 1) -\- ^i{t))/C{t)- Then, the results were grouped according to similarities in both: local 
fields conjugated to the corresponding trade flows, r]ij, Eq. (|13|) . and their year-by-year changes, d^i/^i + d^j/^j, characterizing 
the initial economic conditions influencing trade and the magnitude of the applied perturbation, respectively. Each group of 
connections, V{m,n), was characterized by two integers, m and n. The trade flow, Vij, was classified as belonging to V{m,n) 
if m — 1 < lnr]ij < m and n — 1 < 100 dr]ij < n. (Note that the grouping with respect to m is in fact logarithmic binning 
with respect to the expected trade volume, Eq. (p!3|) . while the parameter m describes linear binning with respect to the sum of 
percentage changes in relative GDPs of trade partners.) The gray circles shown in the figure correspond to geometric averaging 
of the results pertaining to single connections of ITN. The averaging was employed over the predefined groups, V{m,n). The 
black squares correspond to grey points averaged over the horizontal axis. More precisely, the upper left panel a, shows all 
the data such as they are described above. The upper right panel b, presents the same data with the additional condition: it 
illustrates the fluctuation-response theorem as applied to trade volumes whose expected share, (fij), Eq. (|13|). in the global 
trade was greater than 10""^. The lower left panel, c, presents those points from the panel a, for which the number of single 
trade volumes contributing to the geometric average was greater than 1000. Finally, in the lower right panel d, the gray circles 
which meet the additional conditions specified in panels b, and c, are shown. 



